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We study the anisotropic lattice QCD for precision computations of heavy-light matrix elements. Our previous 
study in which the lattices are calibrated with a few percent accuracy has already given results comparable to 
the existing calculations. This suggests that even higher precision may be achieved by a more precise calibration 
of anisotropic lattices. We describe our strategy to tune the gauge and quark parameters with accuracies much 
less than 1 % in the quenched approximation. 



1. Introduction 

Recent experimental progress at B factories 
suggests that precise computation of hadronic 
matrix elements in lattice QCD is a key to the 
search for signals of new physics in flavor physics. 
However, HQET or the relativistic approaches 
still suffer from perturbative and/or discretiza- 
tion errors, which are typically of order of 10%. 
We therefore need yet another framework of the 
heavy quark in which one should be able to (i) 
take the continuum limit, (ii) compute the param- 
eters in the action and the operators nonpertur- 
batively, (iii) and compute the matrix elements 
with a modest computational cost. 

As a candidate of framework which fulfills these 
conditions, we investigate the anisotropic lattice 
on which the temporal lattice spacing a T is finer 
than the spatial one a a |H2i;Sj . The anisotropic 
lattice approach evidently satisfies above con- 
ditions (i) and (iii). Our expectation is that 
on anisotropic lattices the mass dependence of 
the parameters becomes so mild that one can 
adopt coefficients determined nonperturbatively 
at massless limit. For precise computations of 
heavy-light matrix elements, we also need to con- 
trol all the systematic errors in the extrapolations 
to the continuum limit. Whether these promises 
will be practically satisfied should be examined 
numerically. 

So far we have investigated the feasibility of 
the approach in the quenched approximation and 



with the tree-level tadpole improvement for the 
0{a) improved Wilson quark action |1I21 '3 4 . As 
will be summarized in Sec. |3| the results have 
been encouraging for further development in this 
direction. Therefore we have started the sec- 
ond stage of the study of the anisotropic lat- 
tice for heavy quarks. In this stage, we perform 
fully nonperturbative 0{a) improvement and aim 
at 0(2%) computations of matrix elements in 
quenched approximation. To this end, we need 
to perform the calibrations with accuracies much 
less than 1 %. We describe our strategy to tune 
the gauge and quark parameters to this level in 
the quenched approximation. 

2. Anisotropic lattice quark action 

Our heavy quark formulation basically follows 
the Fermilab approach \5\ but is formulated on 
the anisotropic lattices |1I6| . The quark action is 
represented as 



Sf = ^i>(x)K{x,y)il>(y), 
K(x,y) = 6 X . V -K T (1 - ^ 4 )Ui(x)S x+ g iy 
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where K a and the spatial and temporal 

hopping parameters, r the spatial Wilson pa- 
rameter and ce and cb the clover coefficients. 
For a given in principle, the four parameters 
7i? = k t /k ct , r, eg and cb should be tuned so that 
the Lorentz invariance holds up to discretization 
errors of 0(a 2 ). We can set r = l/£ without loss 
of generality |1I5| . 

In the first stage of this work, we tuned only 
the bare anisotropy nonperturbatively and ap- 
plied the tree-level tadpole-improvement to eg 
and cb- ce = 1/uo-u 2 ., cb — 1/ u ct- ln the sec- 
ond stage, however, we need to perform the non- 
perturbative tuning of all the three parameters 
7f, ce, and cb, as well as the parameters of the 
operators which appear in the matrix elements. 



3. Summary of the first stage results 

In the first stage of this work, we have obtained 
the following results in the quenched approxima- 
tion. 

One-loop perturbative calculation 0: Renor- 
malization factors of heavy-light bilinears and 
quark rest mass at rriQa a ~ 1, rriQ <C a^ 1 are cal- 
culated in the one-loop perturbation theory. The 
mQa T dependence of the coefficients are well ap- 
proximated with linear form and this means that 
the quark mass dependence can be controlled. 

Numerical simulation in the quenched approxi- 
mation 0: The mass dependent tuning of 7f is 
performed with meson dispersion relation. Quark 
mass dependence is small for m q a T <C 1 and well 
fitted to a linear form in (m q a T ) 2 . 

Test of relativity relation [3]: Heavy- light me- 
son dispersion relation is computed with j F = 
"f F (m q — 0). The relativity relation well holds 
for the region m q a T <C 1 while m q a a ~ 1. 

Application to the decay constant 0]. Around 
the charm quark mass, the heavy-light decay con- 
stant is calculated. The result is consistent with 
previous works. This implies that our program 
successfully works at least for computations at 
0(10%) accuracy. 

These results are encouraging for further devel- 
opment along this direction. 



4. Strategy for the second stage 

The goal of the second stage is to perform the 
quenched computations of heavy-light matrix el- 
ements within a few percent uncertainties. To 
achieve this precision, the calibration must be 
performed to the level of accuracy much less than 
one percent both for the gauge and quark fields. 

4.1. Calibration of gauge field 

In the quenched approximation, the calibration 
of gauge field can be performed independently of 
the quark field. The elaborated work by Klassen 
0, the 0(1%) level calibration for the Wilson 
action, is no longer enough for the present pur- 
poses. For more precise calibration of gauge field, 
we need to measure the static quark potential 
very accurately. For this purpose, we adopt the 
Liischer-Weisz noise reduction technique 

We define the renormalized anisotropy £g 
through the hadronic radius ro measured in the 
spatial and the temporal directions. Since we 
carry out the continuum extrapolation in terms 
of the lattice scale set by ro, the renormalized 
anisotropy is kept fixed during the extrapolation. 
This procedure prevents the systematic uncer- 
tainties due to the anisotropy from remaining in 
the continuum limit. 

Figurenshows a result of static quark potential 
at (3 — 5.8 and 7g = 3.10: the displayed data are 
the force between static quarks separated in the 
spatial (coarse) and temporal (fine) directions. 
The values of ro, which is defined through the 
relation r^F^o) — 1.65, are determined within 
0.2% statistical errors. Fig. |3 exhibits the deter- 
mination of Jq, which satisfies £g(7g) = £ — 4, at 
(3 = 5.75. A linear X 2 fit gives 7g = 3.1399(52), 
which satisfies the required accuracy. The cali- 
bration in wider range of j3 at £ = 4 is in progress. 

4.2. Calibration of quark field 

We need to calibrate the parameters in the ac- 
tion, tf, ce, cb, to the level which enables com- 
putations of matrix elements within a few percent 
accuracy. We also need to perform the nonper- 
turbative renormalization of the operators such 
as the heavy-light axial current. The nonpertur- 
bative renormalization technique [10111112] is one 
of the most powerful methods to perform such a 



3 



0.5 




0.0 1 ' ' ' 1 ' ' ' 1 ' ' 

0.0 2.0 4.0 6.0 

r. 

Figure 1 . Determination of gauge field anisotropy. 
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Figure 2. Calibration of gauge field at (3 — 5.75. 

program. Since our expectation is that the result 
of tuning near the massless limit can be applied 
to the heavy quark region if rriQa T <C f holds, 
the technique can be applied with a little modifi- 
cation in accord with the anisotropic lattice. 

The Schrodinger functional method can be ap- 
plied to the anisotropic lattice in a straightfor- 
ward manner if the fine direction is realized as 
the temporal axis. To examine the feasibility 
of the method on anisotropic lattice, we perform 
the tree level analysis in nonzero background field 
along Ref. Requiring the PC AC relation up 
to 0(a 2 ), the tree level relations £/7^ = 1 and 
= 1 are reproduced. In this setting, eg can 
be tuned with sufficient accuracy, while Cb seems 
not. The calibration of may also be insuffi- 



cient, since it is the parameter to be tuned most 
precisely. 

To tune all the required parameters to a suf- 
ficient level, we perform the calibration of jp, 
ce, cb, and the renormalization coefficients of the 
axial current along the following steps. (1) Tun- 
ing of Ce by Schrodinger functional method (and 
7f if sufficient accuracy is accessible). (2) Cal- 
ibration of 7f and cb by requiring the physical 
isotropy conditions for mps and my in the coarse 
and fine directions on lattices with T,L ~ 2 fm. 
(3) Determination of k c and the renormalization 
coefficients of the axial current by Schrodinger 
functional method. (4) Finally, several checks are 
necessary. We need to verify that the systematic 
errors are under control by calculating the hadron 
spectra and the dispersion relations and by taking 
the continuum limit. It is also necessary to verify 
that the tuned parameters in the massless limit 
is also available in the heavy quark mass region. 

The numerical simulation along this program 
is in progress. 

REFERENCES 

1. J. Harada et al, Phys. Rev. D 64 (2001) 
074501. 

2. H. Matsufuru, T. Onogi and T. Umeda, Phys. 
Rev. D 64 (2001) 114503. 

3. J. Harada, H. Matsufuru, T. Onogi and 
A. Sugita, Phys. Re v. D 6 6 (2002) 014509. 

4. H. Matsufuru et al, hep-lat/0209090 

5. A. X. El-Khadra, A. S. Kronfeld and 
P. B. Mackenzie, Phys. Rev. D 55 (1997) 3933. 

6. T. Umeda et al, Int. J. Mod. Phys. A 16 
(2001) 2215. 

7. T. R. Klassen, Nucl. Phys. B 533 (1998) 557. 

8. M. Liischer and P. Weisz, JHEP 0109 (2001) 
010. 

9. R. Sommer, Nucl. Phys. B411 (1994) 839. 

10. M. Liischer et al, Nucl. Phys. B478 (1996) 
365. 

11. M. Liischer and P. Weisz, Nucl. Phys. B479 
(1996) 429. 

12. M. Liischer et al, Nucl. Phys. B491 (1997) 
323. 



